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This paper aims to study the multi-level inverter (MLI) fed the widely used 
doubly-fed induction generators (DFIG). The new DFIG generator connected 
to the electrical grid is evaluated during single-pole autoreclosure operation 
to continue the generator support to the system. The effect of using MLI on 
the performance of the DFIG system during the dynamic operation of the 
autoreclosure is studied for the first time and integrating a new MLI during 
the fault and clearing periods. A dynamic arc model is represented depending 
on Kizilcay model. Furthermore, the dynamics performance of the DFIG 
during integration with MLI to improve the performance is studied. The MLI 
is preferred because of processing high power with high resolution of the 
voltage waveforms and lower stresses on switches. The MLI enhances the 
performance of DFIG-based generating system associated with electrical grid 
system in case of single-pole autoreclosure. The proposed cases of studies 
are simulated using Matlab/Simulink to evaluate the proposed system 
performance during the autoreclosure operation for faults in double-circuit 
transmission systems. 
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1. INTRODUCTION 

The Wind energy has obviously become an effective support for the energy sources in electricity 
generation all over the world. The transmission systems are important to connect between generation units 
and loads and accordingly facilitate bulk power transmissions. For example, the wind farms generations are 
delivered to the loads through the transmission systems. Generally, the double-circuit transmission systems 
are used in order to enhance the transmission system reliability. If one of the circuits is faulted, it is isolated 
by the protection systems and accordingly, the power is transmitted through the other circuit. 

Rcently doubly-fed induction generator (DFIG) is one of the resources of wind energy converted to 
electric energy. DFIG is the commonly used pattern-based wind turbine because it has small scale converters 
and its control facility regarding the active and reactive power which is significant regarding the grid 
integration [l]-[4]. Previously, the generators-based wind turbine farms were out of service from the grid 
connection when the terminal voltage less than 80% of the rated value, which is currently not valid because 
of high breakthrough of wind farms. For regulating the wind farms operation, the grid code requirements 
were updated by the power system operation [5]-[7]. The worthiest demand in these regulations is that related 
to the behavior during the fault cases. Through this requirement, the generator immunity should be increased 
during the fault and stayed connected to the grid for a specific period defined by this rule characteristic, 
which is called the fault ride through (FRT) or the low voltage ride through (LYRT) characteristic [8], [9]. 
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There are several fault classifications in power systems such transient/permanent faults, single- 
phase/two-phase/three-phase faults, earth/phase faults, and shunt/series faults. When the fault is permanent, 
the line should be isolated without any reclosing in order to avoid reclosing onto fault. However, when the 
fault is transient, the line should be isolated and then reclosing after a predetermined time in order to avoid 
closing onto fault. Most of the faults in the power transmission systems are transient faults where they 
represent more than 80 %. The main feature of these faults is that they are dielectric breakdown through the 
arc element and accordingly the transient faults can be called arcing faults. 

When the arcing faults are detected using the protection system, the input energy to the arc is 
diminished when the fault current is interrupted. However, there is a time where the arcing medium is de¬ 
ionized where this period is called the secondary arc. The secondary arc period time depends on the 
surrounding medium such as wind. After the arcing medium is de-ionized, the reclosure can safely reclose 
without switching onto fault and accordingly it inserts the circuit onto healthy operation. Generally, the 
reclosing instant depends on detecting the secondary arc extinction instant and the reclosure is adaptively 
reclosing that is called adaptive reclosure [10]-[17]. Integrating the designed DFIG with the autoreclosure 
protection system is evaluated in this paper. However, the operation of the single-pole autoreclosure has not 
integrated with the performance of the grid connected DFIG. 

Using multi-level inverters (MFIs) [18], [19] is preferred compared with the conventional two levels 
H-bridge inverter for many reasons. MFIs are able to process high power at low switching frequency and 
high output voltage resolution (low total harmonic distortion (THD)). Therefore, MFIs have lower switches 
stresses with higher efficiency. One more important advantage of the MFIs is the modularity which enable 
MFIs to operate even if under faulty condition [20], [21]. This is can be done by modifying the control 
algorithm to bypass the faulty section without changing in the MFI structure. 

There are three main MFIs categories (1) neutral point clamped MFI; (2) cascaded dc sources 
inverter, (3) flying capacitor MFI inverter. MFIs can also be divided into two main categories based on the 
dc sources available, symmetrical and asymmetrical MFIs. In symmetrical MFI, all dc sources are equal 
while in asymmetrical MFIs dc sources are not equal. Asymmetrical MFIs are used to generate higher levels 
compare to the symmetrical MFIs for the same number of switches. The main ratios among dc sources are 
1:2 and 1:3, and the 1:3 ratio can develop the highest output voltage resolution and thus lowest THD. 
Recently hybrid MFI is emerged [22]-[24] to reduce the number of MFI’ switches used and thus increasing 
their reliability. Hybrid MFIs have different structure with different dc sources ratio. By looking to the 
literature review, one can notice that MFIs are rarely incorporated with the DFIG system although they have 
many types, categories and advantages. This is indicated that MFI integrated with the DFIG system is open 
for research. 

In this paper, a multi-level inverter is integrated with the DFIG instead of the conventional two 
levels inverters. The performance of multi-level inverter control is evaluated when the wind driven system 
fed DFIG is subjected to transient ground fault. This research reveals that the wind turbine generator 
immunity could be capable of supporting the grid through the single-phase autoreclosure. Then, the DFIG 
during the fault periods in the double-circuit transmission systems and clearing fault is integrated. The 
double-circuit is considered in this paper due to its enhancement in the grid reliability as the faulted circuit is 
interrupted using the protection while the healthy circuit continues the power transmission service. This 
double-circuit interconnection increased the DFIG low voltage ride through capability during the single-pole 
autoreclosure operation. 


2. SIMULATED SYSTEM 

Figure 1(a) shows the simulated power system containing the 380-kV power grid, double-circuit 
transmission line, delta/star transformer and wind farm. For the wind farm installation, the wind speed is 
assumed as 15 m/s. Although the wind has a dynamic time characteristic, it is considered a fixed value during 
the fault as the fault event and clearance has electrical time constant that is very low comparing the wind 
speed time change. The wind farm consisted of six wind generation units where each unit generated 1.5 MW. 
Further parameters are in the Appendix. The fault case occurred at the middle of one of the circuits as 
depicted in the figure. This transmission system interconnected between the grid and the wind farm where the 
grid simulated by a static source behind impedance and power transformer. 

Figure 1(b) shows a schematic diagram of the structure of the DFIG. The three-phase wound 
induction machine is used. Between the rotor and the grid, there are two converters and capacitor where the 
capacitor voltage is the MFI main stage voltage equal to three times of the auxiliary stage voltage (3 Vd C ). 
The rotor side converter is implemented as a conventional one while the grid side converter is replaced by 
multi-level converter. When the power direction is from the rotor to grid, this converter operation is multi- 
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level inverter while it is operated as a rectifier when the power is in the reverse direction through the dc link. 
Building up the multi-level inverter in the Simulink/Matlab program is discussed as follows. 



(a) Grid-interconnected wind farm (b) Schematic diagram of the DFIG 

Figure 1. Simulated benchmark 


The wind farm under study consists of six units connected in parallel, the parameters of one unit 
based DFIG, converters parameters, and controller parameters are illustrated in Tables 3, 4, and 5. The 
description of these parameters was reported in [28] where they were designed for two conventional 
converters of the DFIG and they are exploited for integrating the proposed MLI. 


Table 3. Parameters of DFIG _ _ Table 4. Converters parameters 


Parameter 

Symbol 

Value 

Parameter 

Value 

Power 

P 

1.5 MW 

Grid side converter coupling inductor 

0.3 p.u. 

Stator voltage 

Vs 

575 Vrms 

Grid side converter coupling resistor 

0.003p.u. 

Rotor voltage 

Vr 

1975 Vrms 

DC Bus voltage (3 Vdc) 

700 V 

Stator Frequency 

F s 

60Hz 

DC Bus capacitor 

0.01 F 

Stator resistance 

R s 

0.023 p.u. 

Line filter capacitor (Q=50) (var) 

120e3 

Rotor resistance 

Rr 

0.016 p.u. 



Rotor leakage inductance 

U 

0.16 p.u. 



Stator leakage inductance 

Lis 

0.18 p.u. 

Table 5. Controller parameters 


Mutual inductance 

Lm 

2.9 p.u. 

Parameter 

Value 




DC bus voltage regulator gains [K p Ki] 

[4 400] 




Grid-side converter current regulator gains [K p 

[0.83 5] 




Ki] 





Rotor-side converter current regulator gains 

[0.6 8] 




[Kp Ki] 





Speed regulator gains [K p Ki] 

[1.5 0.6] 




Q and V regulator gains [Ki var Ki volt] 

[1 30] 


2.1. Mathematical Modeling of Interacted DFIG 

The mathematical model of the DFIG is represented as in [2]. The direct and quadrature axis 
(dq-axis) representation of the machine voltages equations and flux linkages are as follows. 
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where co s is the angular stator speed, 0 ^ is the angular rotor speed, R r and R s are the rotor and stator 
resistances, respectively. The stator inductance is L s , rotor inductance L r and the magnetizing inductance L m . 

L =L +L , L = L + L ( 5 ) 

s Is m 5 r lr m V / 

where L/ 5 and L/ r are the stator and rotor leakage inductances, respectively. The stator active and reactive 
power components can be described by; 

p , = faL+VjJ , Q s - {v qs i ds -VjJ ( 6 ) 

2.2. Multi-level Inverters Design 

Figure 2 shows the three-phase configuration of the hybrid multi-level inverter topology that can 
generate nine levels for line-to-line voltage [23], [24]. This inverter is composed of single unit as a main 
stage, 3 auxiliary units, 4 18’ switches and ‘4’ isolated unequal dc voltage sources. The ratio of the dc voltage 
sources of the auxiliary stage to the main stage is chosen 1/3, therefore this inverter provides high output 
voltage levels. Therefore, if the main stage has voltage source 3V dc , then the auxiliary stage voltage source is 
V dc . The main stage is a conventional two-level three-phase six switch inverter. Each auxiliary cell consists 
of the well-known 4 switches H-bridge single phase inverter with a single dc source. Each auxiliary cell can 
generate three levels voltages -Vac, 0, and V dc . The main stage can generate two levels either 0V or 3V dc . 
Therefore, the generated phase voltage to the neutral of the dc source v aN will have five levels, they are 0, 
Vac, 2V dc , 3V dc , and 4V dc . Therefore, the load line-to-line voltages can have 9 levels voltages (4V dC: 3V dc , 
2V dc ,V dc ,0, -V dc , -2V dc , -3V dc , -4V dc ,). Table 1 shows the switching states of this inverter topology. The idea 
behind chosen this inverter is that, only single dc source can be used for the main stage and the other 3 dc 
sources for the auxiliary circuits can be generated using a high frequency transformer (HFT) [25]. 



Table 1. Switching States of phase v aN 
Switches of arm ‘A’ 


v aN 

Sal Sa2 

Sa3 

Sa4 

Sa5 

Sa6 

OVdc 

1 0 

1 

0 

0 

1 

lVdc 

1 0 

0 

1 

0 

1 

2 Vdc 

0 1 

1 

0 

1 

0 

3V dc 

1 0 

0 

1 

1 

0 

4 V dc 

1 0 

0 

1 

1 

0 


Figure 2. Power circuit of the nine-level line-to-line 
hybrid MLI 


2.3. Space Vector Control Algorithm of Three-phase Asymmetrical Cascaded MLI 

The inverter output voltages are generated using space vector modulation control scheme. The ac 
reference voltages are approximated by three nearest switching voltages by the inverter. Therefore, this 
section presents the well-known space vector modulation (SVM) control algorithm [22]-[24]. Each voltage of 
v aN > v bN and v c n has five possible states (0, V dc , 2V dc , 3V dc , 4V dc ) as shown in table 1, thus there are 125 
possible states of this inverter generating 61 voltage vectors shown in Figure 3. These 61 vectors can be 
divided into four groups based on their magnitude. It is observed from Figure 3 and Table 2 that each vector 
group located on single hexagon constituting 4 groups. 
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(a). Group 2: V ± to V 6 
with V dc magnitude 


(b). Group 3: vectors 

Vy tO Via with ^/3V dc 
magnitudes (solid) and 
2V dc (dash) 


(c). Group 4: V 19 to V 36 
with 2.64V dc (solid) and 
3V dc (dash) 


(d). Group 5: V 37 to V 61 
with 2A6V dc (solid) and 
4V dc (dash) 


Figure 3. Voltage vector groups diagram of the five-level 3-phase MLI 


Table 2. Groups of space voltage vectors 


Vector name 

States 

magnitude 

Group 1: Zero vectors V 0 

[000], [111], [222], [333] and [444]. 

OV 

Group 2: Small vectors V 1 A V 6 

Six vectors with 24 states 

v dc 

Group 3: Medium vectors 1 V 7 A V 12 

Six vectors with 18 states 

V3V dc 

Group 3: Medium vectors 2 V 13 -> V 18 

Six vectors with 18 states 

2V dc 

Group 4: Medium vectors 3 V ig -> V 30 

Twelve vectors with 24 states 

2.65V dc 

Group 4: Medium vectors 4 V 31 A V 36 

Six vectors with 12 states 

3V dc 

Group 5: Large vectors 1 V 37 A V 42 

Six vectors with 6 states 

3A6V dc 

Group 5: Large vectors 2 V 43 -> V 54 

Twelve vectors with 12 states 

3.6lV dc 

Group 5: Large vectors 3 V 55 -> V 61 

Six vectors with 6 states 

4V dr 


In this S VM method, the normalized space vector for every switching state states for the three-phase 
voltages v aN ,v bN and v cN with respect to the neutral of the dc bus voltage is defined by; 

/ n / 12E 1211 \ 

Vsvn = [—) * y v aN + e 3 v bN + e 3 v cN J = V svd + j V SVq (7) 

The magnitude of this normalized vector and its angle can be governed by; 

IlKsvwll = JvJvd + V svq and P = tan ~ l Q ( four quadrant angle) (8) 

The magnitude of the normalized space vector defines the group of vectors that the space vector lays. These 
steps of defining the group of vectors are provided below. Once the group of vectors is defined then the angle 
of the normalized space vector is easily used to locate the normalized vector inside the triangle of the three 
nearest vectors (F 0 , V 1 , V 2 ). Finally, the switching times (T 0 ,7\ , T 2 ) of the three nearest vectors are 
accurately defined from; 


T 0 *V 0 + T 1 *V 1 + T 2 *V 2 - T * Vsvn (9) 

where T = — and is the sampling interval. From table 2, the following steps summarize how to implement 

fs 

SVM control algorithm. 

Step 1: If IlF^wll <0.5 then V SVN lays in Group 1 (zero voltage vector). 

Step 2: If 0.5 < HVs^ll ^ 1.366 then V SVN lays in Group 2. 

Step 3: If 1.366 < ||V 5 ra || < 2.825 then V SVN lays in Group 3. 

Step 3: If 2.825 < HV^^H <3.23 then V SVN lays in Group 4. 

Step 5: If ||Fsy^v || > 3.23 then V SVN lays in Group 5. 

2.4. Arcing Fault Model for Single-pole Autoreclosure 

The single-pole autoreclosure has been incorporated in different power transmission systems as it 
enhanced the system stability. Several challenges are the motivation of research continuity such as 
discrimination between transient and permanent faults, secondary arc extinction detection to avoid reclosing 
onto fault, and accurately estimating the reclosing instant to provide the best reclosing system performance 
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[ 10]-[ 17]. However, such studies and others did not evaluate the effect of single-pole autoreclosure operation 
on the DFIG performance. Also, incorporating the autoreclosure in the double-circuit transmission systems is 
not evaluated for the wind driven DFIG. 

Reference to the simulated system shown Figure 1 and as this fault is considered as a transient fault 
type, the arcing element between the conductor and ground is implemented using the dynamic arc by [26]; 


dg_ 

dt 


if H\ 

T\(U 0 +r\i\)l & 


( 10 ) 


Where g is the time dependent arc conductivity, t is the time, r is the arc time constant, i is the arc 
current, / is the arc length, r is the arc resistance per unit length, and U Q is the arc voltage per unit length. In 
order to implement the arcing fault, the dynamic arc model (10) is solved for each instantaneous arc current i. 
This solution is done with defining the arc parameters concerning the Kizilcay arc modeling that the arc time 
constant is 1.3 ms, arc voltage per unit length is 12 V/cm, arc resistance per unit length is 1.3 mQ/cm, and arc 
length is 350 cm as defined in [10], [27], [28]. This arc implementation is suitable only for the primary arc 
period that has characteristics with high arcing current. However, for the secondary arc period, the arcing 
current is extremely reduced due the single-pole breaker opening at both terminals of the faulted phase. 
Although the same differential equation is used, the arc parameters are changed. However, the secondary arc 
period does not affect the dynamic performance of the DFIG and the DFIG sees the system as it is in normal 
operation with missing single-phase in the double circuit transmission system. Accordingly, this secondary 
arcing period is not represented in the simulation benchmark. Only, the period of 0.24 sec observed in the 
simulation in [10], [26], [27] is considered as the time sufficient for extinguishing the secondary arc period. 
Then, the system can restore reclosing the isolated circuit to the normal operation. The tower resistance is 
considered 12 Q that is added in series with the arcing model. 


3. RESULTS AND ANALYSIS 

Integrating the models of the DFIG considering the data given in the Appendix, MFI, SVC, dynamic 
arc model, and double transmission line is accomplished using the Matlab/Simulink program. Accordingly, 
the simulation-based evaluation of incorporating the MFI within the DFIG instead of one of the conventional 
converters is evaluated with the single-pole autoreclosure as follows. 

3.1. Arcing Fault Evaluation 

First, the dynamic arc characteristics are evaluated as shown in Figure 4. Referring to the system 
shown in Figure 1(a), the fault occurred at the middle of a circuit in the double-circuit transmission line 
where the fault instant is 1 ms. Accordingly, the fault current is suddenly increased as shown in Figure 4(a) 
while the fault is reduced at 1.202 s when a single-pole of the breaker CBa 2 operates to isolate the fault from 
the grid side. Then, the fault is fully isolated when the single-pole of breaker CBb 2 operates at 1.304 s at the 
DFIG side. Accordingly, there are two different arcing fault periods as depicted by the dynamic arc 
characteristics in Figure 4(b). This performance confirms the bilateral interaction between the arc model and 
the power system. 


i ■ 10< 



(a) Fault current 
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(b) Dynamic arcing characteristics 

(c) 

Figure 4. Arcing fault case at the middle of a circuit in double-circuit transmission system 


3.2. Evaluation of DFIG integrated MLI 

As aforementioned, the autoreclosure scenario is that the fault occurred at 1 s, a single-pole of 
breakers CBa 2 and CBb 2 respectively opened at 1.202 s and 1.304 s to isolate the faulted phase, then the 
secondary arcing fault period was extinguished at 1.54 s, and finally the opened single-poles of the breakers 
CBa 2 and CBb 2 respectively closed at 1.55 s and 1.65 s in order to restore the system normal operation. 
Accordingly, the interaction of the DFIG with the power system transient due to autoreclosure operation is 
evaluated as shown in Figure 5. In Figure 5(a), the fault from 1 s to 1.304 s period produced a slight effect 
although the fault current is high. This is because the transformer delta/star connection opens the zero- 
sequence loop and accordingly the fault current distribution is changed to be as phase-to-phase fault with 
lower current amplitudes as depicted in Figure 5(b). Figure 5(c) shows the DFIG power where there is a 
power swing started just after the fault clearance at 1.304 s. This power swing affects the current flow in 
which there is an oscillation in the current envelope as depicted in Figure 5(b). After the system restoration of 
the normal operation at 1.65, the power and currents are directed to the normal operation. In the voltage 
waveforms in Figure 5(a), there are transient at instant 1.55 s that was due to closing the breaker CBa 2 as the 
first step to connect the isolated circuit in the double circuit transmission system. These travelling waves can 
be reduced by controlling the breaker switching time. Figure 5(d) shows the dc voltage across the capacitor in 
between the two converters. The voltage behavior has slight transients due to the autoreclosure operation. 

The performance of the MLI under the previous autoreclosure scenario is observed via the line-to- 
line voltage depicted in Figure 5(e). It could be noticed here that the ratio of the dc voltages between the 
main stage and auxiliary stages are done using the HFT concept [25]. The line-to-line voltage is obviously 
nine levels as stated before. The MLI is successfully integrated with the DFIG. That behavior shown in 
Figure 5 confirms autoreclosure ride-through capability of the DFIG using the nine levels hybrid MLI. 


u 



i 14 14 i j 14 ii U IT 

Tims PmgI 


(a) 3-phase voltages at PCC 



1 1.1 1J O 1A IS iJt 1.T 

TVMtMCl 


(b) 3-phase currents at PCC 
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(c) Generator active power 




(d) 3V dc voltage 



i 1.1 ' :■ U U (J U i.I 

T 


(e) MLI output ac voltage 

Figure 5. (a) 3-phase voltages at PCC, (b) 3-phase currents at PCC, (c) Generator active power, (d) 3Vd C 

voltage, (e) MLI output ac voltage 


4. CONCLUSIONS 

The hybrid multi-level inverter has been designed and incorporated in the DFIG instead of the 
conventional converter at the grid side. This DFIG integrated MLI has been done for the first time. The 
designed MLI has been nine levels that facilitated the DFIG performance with the single-pole autoreclosure 
operation. The arcing fault was modeled using the Kizilcay dynamic arc model. This fault was happened in 
the double-circuit transmission systems. The provided results confirmed that the DFIG integrated MLI 
enhanced the single-pole autoreclosure ride through capability. 
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